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Abstract. In this paper we continue our study of the Kelvin-Helmholtz (KH) instability in relativistic planar jets 
following the long-term evolution of the numerical simulations which were introduced in Paper I. The models have 
been classified into four classes (I to IV) with regard to their evolution in the nonlinear phase, characterized by 
the process of jet/ambient mixing and momentum transfer. Models undergoing qualitatively different non-linear 
^ . evolution are clearly grouped in well-separated regions in a jet Lorentz factor/jet-to-ambient enthalpy diagram. 

' Jets with a low Lorentz factor and small enthalpy ratio are disrupted by a strong shock after saturation. Those with 

, a large Lorentz factor and enthalpy ratio are unstable although the process of mixing and momentum exchange 

1/^ ■ proceeds to a longer time scale due to a steady conversion of kinetic to internal energy in the jet. In these cases, 

' the high value of the initial Lorentz seems to prevent transversal velocity from growing far enough to generate the 

, strong shock that breaks the slower jets. Finally, jets with either high Lorentz factors and small enthalpy ratios 

JnI ■ or low Lorentz factors and large enthalpy ratios appear as the most stable. 

In the long term, all the models develop a distinct transversal structure (shear/transition layers) as a consequence 
. of KH perturbation growth. The properties of these shear layers are analyzed in connection with the parameters 

CL(' of the original jet models. 
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^ 1. Introduction has been determined for the full set of simulated models, 

•'— j ^ , . 1 1 1 r on average. The agreement between the linear stability 

k> , In the previous paper (Ferucho et al. 2004), hereafter as , . , • i ■ i r t^tt ■ ^ i -i-^ ■ 

I^N( ^ 1 r 1 1 • 1 • analysis and numerical simulations of KH instability m 

Faper 1, we nave performed several simulations and char- , , , . , ^ i ■ i t ^ 

. ,, rr ri-..i • the linear range has been achieved for a very high radial 

'■^ acterized the eftects ot relativistic dynamics and thermo- , r / r, / r, ■ ,^ • . i- \ i • i 

■ " " , .... , PT^TT- , -i- • • 1 resolution of 400 zones/it, (Rj is the jet radius), which 
dynamics m the development of KH instabilities m planar 



relativistic jets. We performed a linear stability analysis 
and numerical simulations for the most unstable first re- 
flection modes in the temporal approach, for three differ- 



appears to be especially relevant for hot jets. The veri- 
fied high accuracy of the simulations made it possible to 
extend the analysis up to nonlinear phases of the evolu- 

. , ^ _ , , r tion of KH instability. We identified several phases in the 

ent values of the Lorentz factor 7 (5, 10 and 20) and a few , ^. r ^^ ^^ i i t ^ ^- i • • 

. / ni - evolution of all the models: linear, saturation and mixing 

difterent values of specific internal energy ot the ict matter , mi r ; i ■ r ,^ 

,^ „ _ - „n „ % ■ 1 1- 1 1 • \ phases, ihe further analysis of the mixing process, oper- 

(from 0.08 to 60.0c ; c is the light speed m vacuum). ^- ^ ■ ^^ ^ ^ i -n i r i ■ 

^ , , ,. atmg during the long-term evolution, will be performed m 

In Paper I we focused on the linear and postlinear 

. 1 r 1 the present paper. 



regime, especially on the comparison of the results of the 
linear stability analysis and numerical simulations in the 
linear range of KH instability. 

We demonstrated first that, with the appropriate nu- 
merical resolution, a high convergence could be reached 
between the growth rate of perturbed modes in the sim- 
ulations and the results of a linear stability analysis per- 
formed in the vortex sheet approximation. A 20% accuracy 



We have found that in each of the examined cases the 
linear phase always ends when the longitudinal velocity 
perturbation departs from linear growth. Then the longi- 
tudinal velocity saturates at a value close to the speed of 
light. This limitation (inherent to relativistic dynamics) is 
easily noticeable in the jet reference frame. We also noted 
a saturation of the transversal velocity perturbation, mea- 
sured in the jet reference frame, at the level of about 0.5c. 
Send offprint requests to: M. Perucho The transversal velocity saturates later than the longitu- 
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dinal velocity at a moment which is close to the peak value 
of pressure perturbation. 

Therefore we concluded that the relativistic nature of 
the examined flow is responsible for the departure of the 
system from linear evolution, which manifests as the limi- 
tation of velocity components. This behaviour is consistent 
with the one deduced by Hanasz (1995, 1997) with the aid 
of analytical methods. 

Our simulations, performed for the most unstable first 
reflection modes, confirm the general trends resulting from 
the linear stability analysis: the faster (larger Lorentz fac- 
tor) and colder jets have smaller growth rates. As we men- 
tioned in the Introduction of Paper I, Hardee et al. (1998) 
and Rosen et al. (1999) note an exception which occurs 
for the hottest jets. These jets appear to be the most sta- 
ble in their simulations (see also the simulations in Marti 
et al. 1997). They suggest that this behaviour is caused 
by the lack of appropriate perturbations to couple to the 
unstable modes. This could be partially true as fast, hot 
jets do not generate overpressured cocoons that let the jet 
run directly into the nonlinear regime. However, from the 
point of view of our results in Paper I, the high stability 
of hot jets may have been caused by the lack of radial 
resolution that leads to a damping in the perturbation 
growth rates. Finally, one should keep in mind that the 
simulations performed in the aforementioned papers only 
covered about one hundred time units (Rj/c), well inside 
the linear regime of the corresponding models for small 
perturbations. In this paper, the problem of the stability 
of relativistic cold, hot, slow and fast jets is analyzed on 
the basis of long-term simulations extending over the fully 
nonlinear evolution of KH instabilities. 

Finally, we found in Paper I that the structure of the 
pertubation does not change much accross the linear and 
saturation phases, except that the oblique sound waves 
forming the perturbation became steep due to their large 
amplitude. In this paper we show that the similarity of all 
models at the saturation time, found in Paper I, will lead 
to different final states in the course of the fully nonlinear 
evolution. 

The plan of this paper is as follows: In Section |21 we 
recall the parameters of the simulations discussed in Paper 
I and in this paper along with a new complementary set 
of simulations introduced here with the aim of making 
the analysis more general. In Section we describe our 
new results concerning the long-term nonlinear evolution 
of KH modes and we discuss our results and conclude in 
Section El 

2. Numerical simulations 

In this paper we continue the analysis of simulations pre- 
sented in Paper I. The parameters of the simulations are 
listed in Tabled The values of the parameters were chosen 
to be close to those used in some simulations by Hardee 
et al. (1998) and Rosen et al. (1999) and to span a wide 
range in thermodynamical properties as well as jet flow 
Lorentz factors. In all the simulations of Paper I, the den- 



sity in the jet and ambient gases are poj =0.1, poa = 1 
respectively and the adiabatic exponent Tj^a = 4/3. 

Since the internal rest mass density is fixed, there are 
two free parameters characterizing the jet equilibrium: 
Lorentz factor and jet specific internal energy displayed 
in columns 2 and 3 of Table ^ Models whose names 
start with the same letter have the same thermodynami- 
cal properties. Jet (and ambient) specific internal energies 
grow from models A to D. Three different values of the 
jet flow Lorentz factor have been considered for models B, 
C and D. The other dependent parameters are displayed 
in columns 5-11 of Table Note that given our choice 
of poj , the ambient media associated with hotter models 
are also hotter. The next three columns show the longitu- 
dinal wavenumber together with oscillation frequency and 
the growth rate of the most unstable first reflection (body) 
mode. The following three columns display the same quan- 
tities in the the jet reference frame. The next two columns 
show the transversal wavenumbers of linear sound waves 
in jet and ambient medium respectively. The last column 
shows the linear growth rate of KH instability in the jet 
reference frame expressed in dynamical time units, i.e. in 
which time is scaled to Rj/cgj- All other quantities in the 
table are expressed in units of the ambient density, poa, 
the speed of light, c, and the jet radius, Rj. 

In order to extend our conclusions to a wider region 
in the initial parameter space, we have performed a new 
set of simulations, which will be discussed only in some 
selected aspects. The initial data for these new simula- 
tions are compiled and shown in the lower part of Table 
n The external medium in all cases is that of model A05. 
From models F to H, internal energy in the jet is increased 
and rest-mass density decreased in order to keep pressure 
equilibrium, whereas the jet Lorentz factor is kept equal 
to its value in model A05. 

The initial momentum density in the jet decreases 
along the sequence A05, F, G, H. Simulations I and J 
have the same thermodynamical values as models F and 
G, respectively, but have increasing jet Lorentz factors to 
keep the same initial momentum density as model A05. 
Finally, in simulations K and L we exchange the values of 
the Lorentz factor with respect to those in runs I and J. 

The numerical setup for the simulations described in 
this paper is the same as in Paper I. We only simulate half 
of the jet (x > 0) due to the assumed symmetry of per- 
turbations. Reflecting boundary conditions are imposed 
on the symmetry plane of the flow, whereas periodical 
conditions are settled on both upstream and downstream 
boundaries. The applied resolution of 400 x 16 grid zones 
per jet radius is chosen for all simulations listed in Tabled 
following the tests described in the Appendix of Paper I. 
In the Appendix of this paper we present a discussion of 
the influence of different longitudinal resolutions on the 
simulation results during the long-term nonlinear evolu- 
tion. Other details of the numerical setup are described in 
Paper I. 

The steady model is then perturbed according to the 
selected mode (the most unstable first reflection mode), 
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Model 



P 



n 



rmA05 
B05 
COS 
DOS 
BIO 
CIO 
DIO 
B20 
C20 
D20 



5 
5 
5 
5 
10 
10 
10 
20 
20 
20 



0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 
0.1 



0.08 
0.42 
6.14 
60.0 
0.42 
6.14 
60.0 
0.42 
6.14 
60.0 



0.008 
0.042 
0.614 
6.000 
0.042 
0.614 
6.000 
0.042 
0.614 
6.000 



0.18 0.059 

0.35 0.133 

0.55 0.387 

0.57 0.544 

0.35 0.133 

0.55 0.387 

0.57 0.544 

0.35 0.133 

0.55 0.387 

0.57 0.544 



0.0027 
0.014 
0.205 
2.000 
0.014 
0.205 
2.000 
0.014 
0.205 
2.000 



0.11 
0.14 
0.44 
0.87 
0.14 
0.44 
0.87 
0.14 
0.44 
0.87 



0.11 
0.15 
0.51 
0.90 
0.15 
0.51 
0.90 
0.15 
0.51 
0.90 



5.47 
2.83 
1.80 
1.71 
2.88 
1.83 
1.73 
2.89 
1.83 
1.74 



0.30 
0.69 
2.00 
2.63 
0.50 
1.91 
2.00 
0.46 
1.44 
2.00 



0.20 0.026 

0.49 0.055 

1.60 0.114 

2.18 0.132 

0.41 0.031 

1.72 0.055 

1.81 0.063 

0.39 0.014 

1.37 0.027 

1.91 0.029 



1.32 
2.62 
5.73 
7.02 
3.59 
9.77 
9.67 
6.51 
13.89 
18.11 



7.20 
7.32 
9.98 
11.56 
10.28 
17.67 
16.58 
18.76 
25.38 
31.43 



0.13 
0.28 
0.57 
0.66 
0.31 
0.55 
0.63 
0.28 
0.54 
0.58 



7.08 
6.84 
8.17 
9.18 
9.64 
14.72 
13.47 
17.60 
21.24 
25.68 



0.53 
1.08 
1.07 
0.24 
0.94 
1.49 
0.20 
0.90 
1.28 
0.31 



0.73 
0.79 
1.05 
1.15 
0.90 
1.01 
1.10 
0.81 
0.99 
1.01 



F 
G 
H 
I 
J 
K 
L 



5 

5 

5 
11.7 
15.7 
15.7 
11.7 



0.01 
0.001 
0.0001 

0.01 
0.001 

0.01 
0.001 



0.77 
7.65 
76.5 
0.77 
7.65 
0.77 
7.65 



0.008 
0.008 
0.008 
0.008 
0.008 
0.008 
0.008 



0.41 0.058 

0.55 0.058 

0.57 0.058 

0.41 0.058 

0.55 0.058 

0.41 0.058 

0.55 0.058 



0.0026 
0.0026 
0.0026 
0.0026 
0.0026 
0.0026 
0.0026 



0.018 
0.009 
0.008 
0.018 
0.009 
0.018 
0.009 



0.02 
0.01 
0.01 
0.02 
0.01 
0.02 
0.01 



2.38 
1.78 
1.71 
2.42 
1.81 
2.43 
1.81 



0.46 
0.66 
0.66 
0.30 
0.35 
0.30 
0.30 



0.53 0.14 

0.53 0.15 

0.48 0.15 

0.50 0.07 

0.44 0.058 

0.55 0.054 

0.31 0.069 



1.23 
1.87 
1.95 
1.11 
1.70 
1.17 
1.52 



2.83 
3.22 
3.23 
2.66 
3.05 
2.80 
2.72 



0.70 
0.75 
0.75 
0.82 
0.91 
0.85 
0.81 



2.55 
2.62 
2.57 
2.42 
2.53 
2.55 
2.26 



3.72 
4.99 
4.71 
3.52 
4.16 
3.87 
2.93 



1.70 
1.36 
1.31 
1.99 
1.65 
2.06 
1.47 



Table 1. Equilibrium parameters of different simulated jet models along with solutions of the dispersion relation (23) 
in Paper I, corresponding to fastest growing first reflection mode, taken as input parameters for numerical simulations. 
The primes are used to assign wavenumber and complex frequency in the reference frame comoving with jet. The listed 
equilibrium parameters are: 7 - jet Lorentz factor, poj - rest mass density, Ej and Ea - specific internal energies of jet 
and ambient medium, Csj, Csa - the sound speeds in jet and ambient medium, p - pressure, r/ - relativistic density 
and enthalpy contrasts and Mj - the jet Mach number. All the quantities in the table, except the last column, are 
expressed in units of the ambient density, poa, the speed of light, c, and the jet radius, Rj. Parameters for the new set 
of simulations arc shown in the lower part of the table. 



with an absolute value of the pressure perturbation am- 
plitude inside the jet of = 10~^. This means that those 
models with the smallest pressure, like model A, have rel- 
ative perturbations in pressure three orders of magnitude 
larger than those with the highest pressures, D. However 
this difference seems not to affect the linear and postlinear 
evolution (sec footnote in Sect.|31). 

3. Results 

Following the behaviour of simulated models, we found in 
Paper I that the evolution of the perturbations can be di- 
vided into the linear phase, saturation phase and mixing 
phase. This section is devoted to describing the fully non- 
linear evolution of the modes described in Paper I. Our 
description shares many points with the framework devel- 
oped by Bodo et al. (1994) for the case of classical jets. 

In order to illustrate the growth of perturbations and 
determine the duration of the linear and saturation phases 
in our simulations, we plotted in Fig. 1 of Paper I the 
amplitudes of the perturbations of the longitudinal and 
transversal velocities inside the jet and in the jet reference 
frame, together with the pressure oscillation amplitude. 
We also plotted the growth of the imposed eigenmodes 
resulting from the linear stability analysis. Both veloc- 
ity perturbations are transformed from the ambient rest 
frame to the unperturbed jet rest frame using the Lorentz 
transformation rules for velocity components. 

In Paper I we defined the characteristic times tun and 
tsat as the end of the linear and saturation phases respec- 



tively. At the saturation time, the perturbation structure 
is still close to the structure of the initial perturbation, 
except that the oblique sound waves forming the pertur- 
bation became steeper, leading to the formation of shock 
waves. Finally, nearly all the simulations lead to a sharp 
peak of the pressure oscillation amplitude. The times at 
which this peak appears, ipoakj are equal to or slightly 
larger than the saturation times for different models. The 
definition of tiin, <sat and tpeak has been illustrated in Fig. 1 
of paper I. Table [5] collects the times of the linear and sat- 
uration phases in the different models (columns 2-4) along 
with other characteristic times (defined in the table cap- 
tion and in Subsection l3.1|l . 

3.1. Fully nonlinear evolution: jet/ambient mixing 

The beginning of the mixing phase can be detected by 
the spreading of the tracer contours. This can be seen in 
Fig. ^ where the evolution with time of the mean width 
of the layer with tracer values between 0.05 and 0.95 is 
shown. The times at which the mixing phase starts (tmix) 
are shown in Table [3 Consistently with the width of the 
initial shear layer in our simulations (around 0.1 i?j), we 
have defined imix as the time at which the mixing layer 
exceeds a width of 0.1 Rj. 

For models with the same thermodynamical proper- 
ties, those with smaller Lorentz factors start to mix ear- 
lier (see Table|21). Moreover, according to Fig.Q] the mod- 
els can be sampled in two (or perhaps three) categories. 
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Fig. 1. Evolution of the mean width of the jet/ambient mixing layer with time for all the simulations. Different types 
of lines are used for models with different internal energies: Continuous line: model A; dotted line: model B; dashed 
line: model C; dashed-dotted line: model D. Line thickness increases with Lorentz factor (from 5, thinest line, to 20 
thickest one). A value of 5 Rj for the width of the mixing layer (horizontal dashed line) serves to classify the evolution 
of the different models. 



Models A05, BOS, BIO and DIO have wide (> 5i?j) shear 
layers which are still in a process of widening at the end of 
the simulation. The rest of the models have thiner shear 
layers (< 5i?j wide) which are inflating at smaller speeds 
(0.5 - 2 lO^^c, in the case of models COS, DOS and CIO; 
0.5 - 1.2 lO-^c, in the case of models B20, C20, and D20). 
A deeper analysis of the process of widening of the shear 
layers as a function of time shows that all the models un- 
dergo a phase of exponential growth extending from i^^ix 
to soon after ipoak- 

We note that those models developing wider mixing 
layers are those in which the peak in the maxima of the 
pressure perturbation as a function of time, Apeak, reach 
values of the order of 70 — 100, with the exception of model 
B20 that has a thin mixing layer at the end of the simu- 
lation but has Apeak ~ 100, and model DIO, which does 
develop a wide shear layer but for which Apeak remains 
small. We also note that (with the exception of model 
DIO) the models developing wide mixing layers are those 
with smaller internal energies and also relatively smaller 
Lorentz factors. 

There are two basic mechanisms that contribute to the 
process of mixing between ambient and jet materials. The 
first one is the deformation of the jet surface by large 
amplitude waves during the saturation phase. This defor- 
mation favors the transfer of momentum from the jet to 



the ambient medium and, at the same time, the entrain- 
ment of ambient material in the jet. From Table |21 it is 
seen that the process of mixing and momentum exchange 
overlap during the saturation phase. 

The second mechanism of mixing starts during the 
transition to the full non-linear regime and seems to act 
mainly in those models with large Apeak- As we shall see 
below, this large value of Apeak is associated with the gen- 
eration of a shock at the jet/ambient interface at tpcak, 
which appears to be the responsible of the generation of 
wide mixing layers in those models. Figure El shows a se- 
quence of models with the evolution of mixing in two char- 
acteristic cases, BOS and DOS, during the late lapse of the 
saturation phase. The evolution of model BOS is repre- 
sentative of models AOS, BOS and BIO. Models B20, COS, 
CIO, C20 and DIO have evolutions closer to model DOS. 
As it is seen in Fig. |21 in the case of model BOS (left col- 
umn panels), the ambient material carves its way through 
the jet difficulting the advance of the jet material which 
is suddenly stopped. The result is the break-up of the jet. 
In model DOS, (right column panels), matter from the jet 
at the top of the jet crests is ablated by the ambient wind 
forming vortices of mixed material filling the valleys. 

The large amplitude of Apeak reached in models AOS, 
BOS and BIO is clearly associated with a local effect oc- 
curring in the jet /ambient interface (see second panel in 
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Model 


ilin 


^mcx 


^mix 




^pcak 


^pcak 


^fmix 




A05 


180 


335 


335 


380 


380 


100 


380 





B05 


125 


175 


185 


200 


205 


70 


210 


215 


COS 


100 


115 


120 


125 


130 


5 


> 595 


195 


D05 


105 


115 


115 


120 


130 


5 


> 595 


185 


BIO 


235 


375 


335 


380 


385 


100 


445 





CIO 


210 


245 


240 


245 


250 


10 


> 595 




DIO 


180 


220 


215 


225 


225 


10 


350 


345 


B20 


450 


775 


625 


760 


780 


100 


> 1000 




C20 


270 


675 


595 


645 


775 


5 


> 1000 


> 1000 


D20 


350 


465 


450 


480 


500 


10 


> 1000 


> 1000 



Table 2. Times for the different phases in the evolution 
of the perturbed jet models, tn^: end of linear phase (the 
amplitudes of the different quantities are not constant any 
longer) . igat ■ end of saturation phase (the amplitude of the 
transverse speed perturbation reaches its maximum), t^i^'- 
the tracer starts to spread, tpcak'- the peak in the amplitude 
of the pressure perturbation is reached, ifmix^ external ma- 
terial reaches the jet axis, t^cx'- the jet has transferred to 
the ambient a 1% of its initial momentum, t^cq'- longitudi- 
nal momentum in the jet and the ambient reach equiparti- 
tion. Apeak: relative value of pressure oscillation amplitude 
at the peak (see Fig. 1 of Paper I). Note that, as a general 

trend, tun <^ ^mcx ~ ^mix ^ ^sat ^ ^pcak ^ ^fmix ^ ^nieq- 



the left column of Fig. |21 and Fig. that leads to the 
jet disruption. The sequence of events preceding the jet 
disruption includes the formation of oblique shocks at the 
end of the linear/saturation phase, the local effect in the 
jet/ambient interface just mentioned and then a super- 
sonic transversal expansion of the jet that leads to i) a 
planar contact discontinuity (see last panels in the left col- 
umn of Fig.|2J), and ii) the formation of a shock (see below) 
that propagates traversally (see Sect. ?). It appears that 
contrary to the velocity perturbations in the jet reference 
frame (see Fig. 1 of Paper I), the maximum relative am- 
plitudes of pressure perturbation are strongly dependent 
on physical parameters of simulations^. The origin of the 
shock in models with Apeak ~ 100, that enhances the tur- 
bulent mixing of the jet and ambient fluids, can be found 
in the nonlinear evolution of Kelvin-Helmholtz instability 
that leads to significant changes of local flow parameters 
at the end of saturation phase. For instance, the oblique 
shock front resulting from the steepening of sound waves, 
during the linear and saturation phases (see Figs. 3-6 of 
Paper I), crosses the initial shear layer at the interface of 
jet and ambient medium. The formation of such oblique 
shock implies a sudden and local growth of gradients of 
all the dynamical quantities at the jet boundary. The local 

^ In order to see how much this peak in relative pressure 
amplitude was influenced by initial relative amplitude to back- 
ground values, we repeated the simulation corresponding to 
model D05 with the same initial relative amplitude as model 
A05 (i.e., three orders of magnitude larger). Results show that 
there is not a significant difference in the peak values of the 
pressure and that the evolution is basically the same as the one 
of the original model. 



conditions changed by these oblique shocks may become 
favourable for the development of other instabilities like 
those discussed by Urpin (2002). 

The generation of the shock wave at the jet/ambient 
interface is reflected in the evolution of the maxima 
of Mj^± = jj.±Vj^j_/{'yc,^Csa) and 7c,„ being the 

Lorentz factors associated with Vj^± and Cga, respectively) 
representing the transversal Mach number of the jet with 
respect to the unperturbed ambient medium. This quan- 
tity becomes larger than 1 around ipeak in those models 
with Apeak ~ 100 (see Fig.0J) pointing toward a supersonic 
expansion of these jets at the end of the saturation phase. 
The fact that in all our simulations the ambient medium 
surrounding colder models (i.e.. A, B) are also colder (see 
Sect. 12) favors the generation of shocks in the jet/ambient 
interface in these models. On the other hand, in the case of 
models with the highest Lorentz factors (B20, C20, D20) 
the transversal velocity can not grow far enough to gen- 
erate the strong shock which breaks the slower jets. 

In order to extend our conclusions to a wider region 
in the initial parameter space, we have performed supple- 
mentary simulations (namely F, G, H, I, J, K, L) with 
the aim of clarifying the effect of the ambient medium in 
the development of the disruptive shock appearing after 
saturation. In all the cases, the external medium is that 
of model AOS. 

From models F to H, internal energy in the jet is in- 
creased and rest-mass density decreased in order to keep 
pressure equilibrium, whereas the jet Lorentz factor is kept 
constant to its value in model A05. The transversal Mach 
numbers at the peak in these models reach values very sim- 
ilar to that of model A05 (~ 14). In fact, the formation of 
a shock at the end of the saturation phase is observed in 
these models as it is in model A05. 

The initial momentum density in the jet has decreased 
along the sequence AOS, F, G, H. Simulations I and J 
have the same thermodynamical values as models F and 
G, respectively, but have increasing jet Lorentz factors to 
keep the same initial momentum density as model AOS. 
In the case of model I we find the same behavior as in 
previous ones: large transversal Mach number, shock and 
disruption. On the contrary, model J behaves much more 
like model B20, with a value of transversal Mach number 
slightly larger than one, strong expansion and almost no 
mixing. Finally, in order to know to which extent this 
change in behavior was caused by the increase in Lorentz 
factor or in specific internal energy, in simulations K and 
L we cross these values with respect to those in I and 
J. The results from these last simulations show that the 
evolution of model M is very close to that of model I: the 
large value of the transversal Mach number at the peak, 
shock and strong mixing; and that the evolution of model 
K is close to that of model J: a weaker shock, expansion 
and no mixing. 

In the models developing a shock, mixing is associated 
with vorticity generated after the shock formation and, in 
the case of models AOS, BOS, BIO as well as F, G, H, I 
and L matter from the ambient penetrates deep into the 
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Fig. 2. Evolution of the jet particle fraction showing the development of mixing in two representative models. Left 
column (model BOS): the ambient material carves its way through the jet difiiculting the advance of the jet material 

which is suddenly stopped. Right column (model D05): the amount of ambient matter hampering the jet material is 
smaller and matter from the jet at the top of the jet crests is ablated by the ambient wind. 
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Fig. 3. Pressure distribution at the onset of the jet/ambient surface distortion at the end of the saturation phase for 
models A05, B05, BIO and B20. The corresponding times are 355, 190, 370 and 765 Rj/c. In the case of models A05, 
B05 and BIO this distortion leads to the formation of a shock. 




Fig. 4. Time evolution of the maxima of the transversal Mach number of the jet with respect to the unperturbed 
ambient medium, Mj.x- See text for further explanations. Lines are as in Figure^ 
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jet as to reach the jet axis (jet break-up). The times at 
which this happens for the different models are displayed 
in Table [2 Note that in models A05 and B05 the en- 
trainment of ambient matter up to the axis occurs just 
after tp^ak, whereas in model BIO occurs later, probably 
because the shock in this model is weaker. The process 
of mixing can be affected by resolution, as small resolu- 
tions may suppress the development of turbulence. We 
have analyzed this in the Appendix in which we focus on 
the influence of longitudinal resolution. 



3.2. Fully nonlinear evolution: jet/ambient momentum 
transfer 



Let us now analyze the evolution of the longitudinal mo- 
mentum in the jets as a function of time. Figure shows 
the evolution of the total longitudinal momentum in the 
jet for the different models. Jets in models B05, COS, D05 
and DIO (also D20) transfer more than 50% of their initial 
longitudinal momentum to the ambient, whereas models 
AOS, BIO and CIO (also B20) seem to have stopped the 
process of momentum transfer retaining higher fractions 
of their respective initial momenta. Models C20 and D20 
continue the process of momentum exchange at the end of 
the simulations but at a remarkable slower rate (specially 
C20). 

In the case of models F, G, H, I and L, the transfer of 
longitudinal momentum is also very efficient. The reason 
why these models, as well as AOS, BOS and BIO, develop 
wide shear layers and transfer more than S0% of their ini- 
tial momentum to the ambient could be turbulent mixing 
triggered by the shock. In the case of models DIO and 
D20, the processes of mixing and transfer of longitudinal 
momentum proceed at a slower rate pointing to another 
mechanism. The plots of the time evolution of the jet's 
transversal momentum (Fig. for the different models 
give us the answer. Jets disrupted by the shock (as AOS, 
BOS, BIO) have large relative values of transversal momen- 
tum at saturation (> 0.04) that decay very fast afterward 
(AOS is an exception). The peak in transversal momentum 
coincides with the shock formation and the fast lateral ex- 
pansion of the jet at ipcak- Contrarily, models DIO and D20 
have a sustained value of transversal momentum after sat- 
uration which could drive the process of mixing and the 
transport of longitudinal momentum. In these models, the 
originally high internal energy in the jet and the high jet 
Lorentz factor (that allows for a steady conversion of jet 
kinetic energy into internal) make possible the sustained 
values of transversal momentum. Between these two kinds 
of behavior are hot, slow models COS, DOS that do not de- 
velop a shock having, then, thin mixing layers, but trans- 
ferring more than S0% of their longitudinal momentum. 



3.3. Fully nonlinear evolution: classification of the 
models 

Our previous analysis based in the width of the mixing 
layers and the fraction of longitudinal momentum trans- 
fered to the ambient can be used to classify our models: 

— Class I (AOS, BOS, BIO, F, G, H, I, L): develop wide 
shear layers and break up as the result of turbulent 
mixing driven by a shock. 

— Class II (DIO, D20): develop wide shear layers and 
transfer more than S0% of the longitudinal momentum 
to the ambient, as a result of the sustained transversal 
momentum in the jet after saturation. 

— Class III (COS, DOS): have properties intermediate to 
models in classes I and II. 

— Class IV (B20, CIO, C20, J, K): are the most stable. 

Figures I7I10I show the flow structure of the different 
models at the end of the simulations. The following mor- 
phological properties of the members of each class are re- 
markable: 

— Class I: irregular turbulent pattern of the ffow, the 
structure of KH modes still visible on the background 
of the highly evolved mean flow pattern. 

— Class II: a regular pattern of "young" vortices (vis- 
ible in the tracer and specific internal energy distri- 
butions), the structure of KH modes visible. The en- 
hanced transfer of momentum found in the models of 
this class is probably connected to the presence of these 
"young" vortices. 

— Class III: the flow is well mixed, i.e. tracer, internal en- 
ergy and Lorentz factor are smoothed along lines par- 
allel to the jet symmetry plane. Highly evolved vortices 
visible. A fossil of KH modes visible only as pressure 
waves. 

— Class IV: no vortices, no chaotic turbulence, very weak 
mixing, very regular structure of KH modes. 

3.4. Transversal jet structure at late stages of 
evolution 

At the end of our simulations, the models continue with 
the processes of mixing, transfer of momentum and con- 
version of kinetic to internal energy, however they seem to 
experience a kind of averaged quasisteady evolution which 
can be still associated with the evolution of a jet, i.e., a 
coUimated flux of momentum. This jet is always wider, 
slower and colder than the original one and is surrounded 
by a broad shear layer. This Section is devoted to the ex- 
amination of the transition layers in distributions of gas 
density, jet mass fraction and internal energy as well as 
shearing layers in velocity, longitudinal momentum and 
Lorentz factor. 

Let us start by analyzing the overall structure of the 
pressure field for the different models at the end of our sim- 
ulations. Figure^Jshows the transversal, averaged profiles 
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Fig. 5. Evolution of the total longitudinal momentum in the jet as a function of time for all the simulations. Lines 
represent the same models than in Fig. ^ The long-dashed horizontal line serves us to identify those models tranferring 
more than 50% of the initial jet momentum to the ambient. 
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Fig. 6. Evolution of the total transversal momentum in the jet as a function of time for all the simulations. Lines 
represent the same models than in Fig.Ul Left panel: models A05, B05, BIO, COS, D05. Right panel: B20, CIO, C2G, 
DIO, D20. Note the change in both horizontal and vertical scales between the two panels. 



of pressure accross the computational grid. Several com- 
ments are in order. In the case of models A and B, the 
shock formed at the end of the saturation phase is seen 
propagating (at r « 30Rj in the case of model A, and at 
r « 40 — 60i?j in the case of models B) pushed by the over- 
pressure of the post shock state. In the case of models C 
and D, the wave associated with the peak in the pressure 
oscillation amplitude seem to have left the grid (remember 
that in our jet models, hotter jets have also hotter ambi- 
ent media). The most remarkable feature in the pressure 



profile is the depression centered at r « 2Rj in the case 
of models CIO and C20 and at r « 3Rj in the case of 
models DIO, D20. These pressure minima coincides with 
the presence of vortices (clearly seen in models C20 and 
D20 in the corresponding panels of Fig. Inland Fig. I10|) . 
Also remarkable in these plots is the almost total pres- 
sure equilibrium reached by models COS and DOS and the 
overpressure of the jet in model D20. 

As we noted in the previous Section the models evolve 
following four schemes. Jets belonging to the classes I and 
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Fig. 7. Snapshot in the mixing phase of logarithmic maps of pressure, jet mass fraction and specific internal energy 
and non-logarithmic Lorentz factor for model A05. Only the top half of the jet is shown. 
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Fig. 11. Longitudinal averaged profiles of gas pressure 
for all models. Different types of lines are used for models 
with different internal energies: Continuous line: model A; 
dotted line: model B; dashed line: model C; dashcd-dotted 
line: model D. Line thickness increases with Lorentz factor 
(from 5, thinest line, to 20 thickest one). 



II disrupt leading to the dispersion of tracer contours for 
more than five initial jet radii. Model D20 of the II class is 



specific. It has not reached the tracer contour dispersion 
equal to five jet radii, but it clearly follows from Fig. ^ 
that this should happen around t — 1100. The models 
belonging to the classes III and IV do not exhibit the 
dispersion of tracer contours for more than 5 jet radii and 
look different at the end of simulations. 

Fig. El displays, for models BOS and DOS representing 
classes I/II and III/IV respectively, the profiles of relevant 
physical quantities averaged along the jet at the end of 
the simulations. Let us note that different shear (in case 
of velocity related quantities) or transition layers (in case 
of material quantities) can be defined depending on the 
physical variable used. 

In case of model BOS all the material quantities (tracer, 
density and internal energy) exhibit a wide broadening in 
the radial direction. The distribution of tracer extends up 
to r = 15Rj, which means that jet material has been 
spread up to this radius, with a simultaneous entrain- 
ment of the ambient material into the jet interior. The 
latter effect is indicated by the lowering of the maximum 
tracer value from 1 down to 0.4. The fine structure of 
the tracer distribution displays random variations, which 
apparently correspond to the turbulent flow pattern well 
seen in Fig. |S| (upper set of panels). The curve of internal 
energy is very similar to the one corresponding to tracer, 
however variations are seen up to r = 20i?j . The profile of 
density is wider than the profile of the tracer (the density 
is growing up to r ~ 40i?j, which can be explained by 
the heating of external medium, in the jet neighborhood, 
by shocks associated with outgoing large amplitude sound 
waves and by transversal momentum transmitted to the 
ambient medium via sound waves. Finally, the profile of 
the specific internal energy is consistent with the density 
profile and the fact of the jet being in almost pressure 
equilibrium. 
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Fig. 8. Same as Fig. for models B05 (upper; only top half of the model shown), COS (middle) and DOS (lower). 
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Fig. 9. Same as Fig. ((JJ for models BIO (upper; only top half of the model shown), CIO (middle) and DIO (lower; 
only top half of the model shown). 
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Fig. 12. Transversal averaged profiles of relevant physical quantities at the end of simulations B05 and D05. Left 
column: tracer, / (full line), rest mass density, po (dotted line), specific internal energy e (dashed line) and jet internal 
energy density , e (= poef; dash-dot line). Right column: longitudinal velocity, Vz (full line), Lorentz factor normalized 
to its initial value, 7/70 (doted line) and longitudinal momentum normalized to its initial value, S/Sq. The upper plots 
represent the model B05 and the lower plots represent DOS. Note that the values of e are multiple by 10 for model 
B05 and divided by 10 for model DOS. The values of e for the model DOS are divided by 100. 



The dash-dot curve in Fig. ^] (top left panel) repre- 
sents the (internal) energy per unit volume held in jet 
matter. Such a quantitiy, like the mean Lorentz factors in 
both inner jet and shear layer are of special importance as 
they are directly related to the emission properties of the 
model. Internal energy density in jet particles is related 
to the fluid rest frame synchrotron emissivity, whereas the 
fluid Lorentz factor governs the Doppler boosting of the 
emitted radiation. As seen in the top right panel of Fig. 1121 
the final mean profile of velocity is similar in shape to 
the profile of internal energy, despite the fact that it is 
smoother. Similarly to internal energy, the longitudinal 
velocity variations extend up to r = 20Rj. This can be 
understood in terms of large amplitude, nonlinear sound 
waves, which contribute to the transport of internal and 
kinetic energies in the direction perpendicular to the jet 
axis. The profiles of Lorentz factor and longitudinal mo- 
mentum are significantly narrower. Therefore in case of 
models similar to BOS only the the most internal part, 
up to r < 8Rj of the wide sheared jet, will be Doppler 
boosted, even though the jet material quantities extend 
behind r ~ l5Rj. 

A similar discussion can be performed for the model 
DOS representing the other group of jets, which form a 
shear layer without experiencing the phase of rapid dis- 
ruption. As seen in the bottom panels of Fig. the jet 
of model DOS preserves sharp boundaries between their 
interior and the ambient medium, although both media 
are modified by the dynamical evolution at earlier phases. 
The sharp boundary (transition layer) at r ~ 3.2Rj is ap- 
parent in the plots of all material quantities, i.e. tracer, 
density, internal energy. The thickness of the transition 
layer for all the quantities is comparable to one initial jet 
radius, 10 times narrower than in case of BOS. We note, 
however a smooth change of ambient gas density in the 
range of r ~ 3.2Rj ^ 12Rj. 



It is apparent also that a narrow core of almost un- 
mixed (/ = 1) jet material remains at the center in the 
currently discussed case. The radius of the core is about 
one half of the original jet radius. The core sticks out 
from a partially mixed, relatively uniform sheath and is 
well seen in the plots of tracer and internal energy for that 
model, however it disappears when increasing the resolu- 
tion in the longitudinal direction, as seen in the Appendix. 

Concerning the dynamical quantities, we note that 
there is no sharp jump in the profiles of longitudinal ve- 
locity, Lorentz factor and longitudinal momentum and the 
central core does not appear in profiles of these quantities. 
Significant longitudinal velocities extend up to r ~ 12-Rj 
as in case of density, in contrast to tracer and internal 
energy. As noted previously, the averaged pressure distri- 
bution for model DOS is practically uniform in the whole 
presented range of the transversal coordinate. Therefore 
as in case BOS we can conclude that the variations of den- 
sity in the ambient medium are due to the heat deposited 
by nonlinear sound waves. On the other hand the widths 
of the profiles of the Lorentz factor and longitudinal mo- 
mentum are comparable to those of jet mass fraction and 
specific internal energy. Then the emission of the whole 
jet volume will be Doppler boosted. 

Models BOS and DOS were considered as representative 
cases of models developing shear layers wider (group 1; 
classes I and II) and narrower (group 2; classes III and IV), 
respectively, than 5Rj. Now the question is up to which 
extent the characteristics of the shearing flow of these two 
models are common to the models in the corresponding 
groups. We note that given the large differences between 
the initial parameters of models in classes I and II, on 
one hand, and III and IV, on the other, we do not expect 
a perfect match among the properties of the transversal 
structure in models within the same group. For example, 
whereas models in class I develop wide shear layers due 
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to the action of a strong shock formed at the end of the 
saturation phase, models in class II develop shear layers 
through a continuous injection of transversal momentum 
and the generation of large vortices at the jet ambient 
interface. 

We now investigate relations among the following av- 
eraged quantities in the whole set of models at the final 
state: the dispersion of tracer contours, the typical widths 
of profiles of density, internal energy density in jet matter 
{re), velocity, Lorentz factor (vry) and longitudinal mo- 
mentum {rs) and the peak values of Lorentz factor and 
the longitudinal momentum (iSmax)- In all cases the peak 
values were measured directly, whereas the typical width 
of the profiles were taken as their width at the mean value 
between the maximum and minimum ones. We find that 
for all models of group 1, > ARj, rs > 3.5Rj, r^ > 7Rj 
and S'/5o < 0.2. In case of all models of group 2,rj < 4i?j, 
rs < 3.5i?j, Te < 7Rj and S'max/S'o > 0.2. 

4. Discusion and conclusions 

We have studied the non-linear evolution of the relativis- 
tic planar jet models considered in Paper I. The initial 
conditions considered cover three different values of the 
jet Lorentz factor (5, 10, 20) and a few different values of 
the jet specific internal energy (from 0.08c^ to 60. Oc^). 
The models have been classified into four classes (I to 
IV) with regard to their evolution in the nonlinear phase, 
characterized by the process of mixing and momentum 
transfer. Cold, slow jets (Class I) develop a strong shock 
in the jet/ambient interface at the end of the saturation 
phase leading to the development of wide, mixed, shear 
layers. Hot fast models (Class II) develop wide shear lay- 
ers formed by distinct vortices and transfer more than 50% 
of the longitudinal momentum to the ambient medium. In 
models within this class, the high Lorentz factor in the 
original jet and its high internal energy act as a source of 
transversal momentum that drives the process of mixing 
and momentum transfer. Between these classes we find 
hot, slow models (Class III) that have intermediate prop- 
erties. Finally we have found warm and fast models (Class 
IV) as the most stable. Whether a jet is going to develop 
a strong shock and be suddenly disrupted seems to be en- 
coded in the peak of the pressure oscillation amplitude at 
the end of the saturation phase and the related transversal 
Mach number. 

The above picture is clarifying but is subject to the 
limitations of our choice of initial parameters that was 
restricted to values with poj =0.1 (see Paper I). This re- 
striction together with the initial pressure equilibrium lead 
to a constant jet-to-ambient ratio of specific internal ener- 
gies for all the models, i.e., hotter jets are surrounded by 
hotter ambient media. In order to extend our conclusions 
to a wider region in the initial parameter space, we have 
performed a supplementary set of simulations (F-L) with 
the aim of disentangling the effect of the ambient medium 
in the development of the disruptive shock appearing after 
saturation. Thus, hot, tenuous, slow/moderately fast jets 



(F, G, H, I, L) behave like cold, dense ones in a cold envi- 
ronment (A05, B05, BIO). However, if these hot, tenuous 
jets are faster (J, K), they behave as warm, fast models 
(e.g., CIO, C20, B20). The fact that the initial Lorentz 
factor is high seems to prevent the transversal velocity 
from growing enough to generate the strong shock which 
breaks the slower jets. 

Models undergoing qualitatively different non-linear 
evolution are clearly grouped in well-separated regions in 
a jet Lorentz factor /jet-to-ambient enthalpy diagram (see 
Fig. USI). Models in the lower, left corner (low Lorentz 
factor and small enthalpy ratio) are those disrupted by a 
strong shock after saturation. Those models in the upper, 
left corner (small Lorentz factor and hot) represent a rela- 
tively stable region. Those in the upper right corner (large 
Lorentz factor and enthalpy ratio) are unstable although 
the process of mixing and momentum exchange proceeds 
on a longer time scale due to a steady conversion of kinetic 
to internal energy in the jet. Finally, those in the lower, 
right region (cold/warm, tenuous, fast) are stable in the 
nonlinear regime. 

Our results differ from those of Marti et al. (1997), 
Hardee et al. (1998) and Rosen et al. (1999) who found 
fast, hot jets as the more stable. The explanation given 
by Hardee et al. (1998) invoking the lack of appropriate 
perturbations to couple to the unstable modes could be 
partially true as fast, hot jets do not generate overpres- 
sured cocoons that let the jet run directly into the nonlin- 
ear regime. However, as pointed out in Paper I, the high 
stability of hot jets may have been caused by the lack of 
radial resolution, that leads to a damping in the pertur- 
bation growth rates. Finally, the simulations performed in 
the aforementioned papers only covered about one hun- 
dred time units, well inside the linear regime of the cor- 
responding models for small perturbations. In this paper, 
the problem of the stability of relativistic jets is analyzed 
on the basis of long-term simulations that extend over the 
fully nonlinear evolution of KH instabilities. 

At the end of our simulations, the models continue 
with the processes of mixing, transfer of momentum and 
conversion of kinetic to internal energy, however they seem 
to experience a kind of averaged quasisteady evolution 
which can be still associated with the evolution of a jet, 
i.e., a coUimated fiux of momentum. This jet is always 
wider, slower and colder than the original one and is sur- 
rounded by a distinct shear layer. Hence transversal jet 
structure naturally appears as a consequence of KH per- 
turbation growth. The widths of these shear (in case of 
velocity related quantities) or transition layers (in case of 
material quantities) depend on the specific parameters of 
the original jet model as well as the physical variable con- 
sidered. However, models in classes HI and IV develop thin 
shear layers, whereas the shear layers of models in classes 
I and II are wider. The possible connection of these results 
to the origin of the FRI/FRII morphological dichotomy of 
jets in extended radiosources will be the subject of further 
research. Extensions of the present study to models with a 
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Fig. 13. Square root of the jet-to-ambient enthalpy ra- 
tio (see Paper I for definitions) versus jet Lorentz factor. 
Symbols represent different non-linear behaviors: crosses 
stand for shock disrupted jets (cold, slow jets, along with 
tenuous, hot, moderately fast or slow ones); diamonds for 
unstable, hot, fast jets; triangles for relatively stable hot, 
slow, and squares for stable, warm, fast, along with hot, 
tenuous, faster jets. 

superposition of perturbations, cylindrical symmetry and 
three dimensions are currently underway. 
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Appendix A: Influence of numerical resolution in 
the nonlinear evolution 

In order to be aware of the limitations of the resolution 
used in the results for the non-linear regime, we repeated 
model COS (C16 here) with the same transversal resolution 
(400 cells/i?j) and changing the longitudinal resolution. 
First we double it (32 cells/i?j, C32) and then multiply 
it by four (64 cells/i?j,C64), and finally, in order to have 
similar resolutions in both transversal and longitudinal 
directions, we performed this simulation using 256x128 
ceUs/i?j (C128). Models C16, C32 and C64 were evolved 
up to a time larger than 600 i?j/c, whereas model C128 
was stopped a.i t — 375i?j/c. 

Table TKa\ displavs the data corresponding to Table |21 
for models C16, C32, C64, C128. Differences in the du- 
ration of the phases are apparent but not significant. 
Regarding the linear regime, perturbations in models C32 
and C64 grow closer to linear predictions (growth rate 
0.093 c/i?j) in both cases; to be compared with the an- 
alytical value 0.114 c/i?j) rather than C16 (0.085 c/i?j). 



C128 has a slower growth rate (0.073 c/i?j) due to its 
smaller transversal resolution (see Appendix of Paper I). 
Values of pressure perturbation at the peak range from 3 
to 8, increasing generally from smaller to larger longitu- 
dinal resolution. Moreover, transversal relativistic Mach 
numbers are also increasing with resolution (from a value 
of about unity to two), so that we observe a stronger, al- 
though still weak, shock in C64 and C128 than in C16 or 
C32. This is maybe due either to the instability giving rise 
to the shock (see Section rrT|l being better captured with 
increasing resolution or as the result of a smaller numerical 
viscosity. 

Figure lXrl shows the time evolution of the mean width 
of the jet/ambient mixing layer and the total longitudi- 
nal momentum in the jet for model C05 as a function of 
resolution. No noticeable differences are found in the evo- 
lution of the different numerical simulations within the 
linear phase (up to t « 100i?j/c). However, there is a 
clear tendency to develop wider mixing layers and enhance 
momentum transfer in those simulations with higher nu- 
merical resolutions as a result of the reduction of numer- 
ical viscosity. In the case of model C128 the processes of 
jet /ambient mixing and momentum exchange are further 
enhanced by the ratio of longitudinal to transversal resolu- 
tion close to unity which favors the generation of vortices 
in the jet/ambient boundary. The enhancement of mixing 
with numerical resolution and the generation of vortices in 
model C128 is seen in the sequence of Panels for different 
simulations (see published paper) show that models C16 
and C32 are very similar, whereas models C64 and C128 
are totally mixed (see the maxima of the tracer values in 
the scales) and colder (as a result of the enhanced mixing 
with the cold ambient medium). 

Finally, mean transversal profiles (appearing in the 
published paper) of relevant physical quantities show that 
the thin, hot core in model C16 disappears in model C128 
due to the enhanced mixing. In model C32 mixing down to 
the axis occurs after a long process. Models C64 and C128 
develop weak shocks after saturation and suffer sudden 
mixing, which goes on in time, cooling down the remain- 
ing relativistic flow. Transition layers in rest mass density, 
jet mass fraction and specific internal energy are wider 
in model CI 28 due to the enhanced mixing. Longitudinal 
velocity and momentum and Lorentz factor profiles are 
more similar in all the cases although the relativistic core 
is thinner in model CI 28. 
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Fig. A.l. Same plots asn(left) and|31(right) for models C16 (solid line), C32 (dotted line), C64 (dashed line) and 
C128 (dash-dotted line). 
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85 
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105 


3.5 
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150 


C64 


80 
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110 


115 


6.5 


125 


150 


C128 


100 


125 
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125 


8.0 


135 


160 



Table A.l. Times for different phases in the evolution of simulations C16, C32, C64 and C128. See Table |21 for the 
meaning of the times at the table entries. 
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